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We search for non-trivial relativistic solutions of the hydrodynamic equations with quasi-inertial 
flows such as in the Bjorken-like models. The problem is analyzed in general and the known results 
are reproduced by a method proposed. A new class of 3D anisotropic analytic solutions with quasi- 
inertial property is found. An ellipsoidal generalization of the spherically symmetric Hubble flow 
with constant pressure is proposed as a particular case. The relativistic expansion of finite systems 
into vacuum is also described within this class. A region of applicability and possible utilization of 
the new solutions for processes of A+A collisions is discussed. 



I. INTRODUCTION 



The equations of the relativistic hydrodynamics have 
highly nonlinear nature and, therefore, only a few ana- 
lytical solutions are known until now. For the first time 
one-dimensional, or (1+1), analytical solution for Lan- 
dau initial conditions - hot pion gas in Lorentz contracted 
thin disk p| , has been developed by Khalatnikov . The 
equation of state (EoS) was chosen as ultrarelativistic 
one: p = c\e,c^ = 1/3. It is noteworthy that according 
to that solution the longitudinal flows developed to the 
end of hydrodynamic expansion, at freeze-out, are quasi- 
inertial: v w XL/t, Much later, in the papers [3J for 
the same EoS have been found the infinite (1+1) boost- 
invariant solution, for finite systems the similar approach 
was developed in Ref.Q. The property of quasi-inertia 
preserves in these solutions during the whole stage of the 
evolution. Bjorken |f| utilized these solution as the ba- 
sis of the hydrodynamic model of ultra-relativistic A+A 
collisions. 

The spherically symmetric variant of a such kind of 
flows with the Hubble velocity distribution, v = r/t, has 
been considered in Ref . ffj . Some generalization of these 
results was proposed in a case of the Hubble flow for EoS 
of massive gas with conserved particle number in Ref. 7] 
and for the cylindrically symmetric boost-invariant ex- 
pansion with a constant pressure in |8(. 

All these solutions were used for an analysis of ultra- 
relativistic heavy ion collisions. Since the longitudinal 
boost-invariance in a fairly wide rapidity region is not ob- 
served even at RHIC, as it was expected, the Hubblc-likc 
models are also used now for a description of the experi- 
mental data J9|. It is naturally, however, that, unlike to 
the Hubble type flows, the velocity gradients should be 
different in different directions since there is an initial 
asymmetry between longitudinal and transverse direc- 
tions in central A+A collisions and, in addition, between 
in-plane and off-plane transverse ones in non-central col- 
lisions. In this letter we make a general analysis of the 
hydrodynamic equations for the quasi-inertial flows aim- 
ing to find a new class of analytical solutions with 3D 
asymmetric relativistic flows. 



II. GENERAL ANALYSIS 

Let us start from the equations of relativistic hydro- 
dynamics: 

d v T^ = 0, (1) 

where the energy-momentum tensor corresponds to a per- 
fect fluid: 

T» v = {e+ p)uV -p-gV (2) 

We can attempt to find a particular class of solutions and 
therefore have to make some simplifications of QJ. We 
do not fix EoS at this stage. 

• Let us put the condition of quasi-inertiality 

u v d v v? = (3) 

which means that flow is accelerationless in the rest 
systems of each fluid element; this property holds 
for the known Bjorken (boost-invariant) and Hub- 
ble flows. 

Then, we find that u^[{e+p)d v u u + u v d u e] + [<u/ i u v d v p- 
d^p] = 0. Contracting this equation with u M we get 

(e+p)d v u v + u"d v e = Q. (4) 

Obviously, the remaining equation to satisfy is: 

u»u v dvp -d"p = (5) 

The task is to find solution of the system © , an d 
. As one can see, the number of equations exceeds the 
number of independent variables. So, the equations must 
be self-consistent in order to have nontrivial solutions. 

We see that Eq.(0} can be rewritten in the form: 

u^d ll s = -F{e){d v u v ) (6) 
where F(e) — e + p, and Eq.JSJ) as the following: 

p'{e){u»u v d v s-d»e) = 0, (7) 
supposing EoS in the form p — p(e). If p'(e) + 1 

u"F(e)(d l/ u u ) + d' 1 e = 0. (8) 
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Normally F(e) ^ 0, and we can divide the last equation 
by F(e) and introduce the function $(e) by the definition 

i+fey = so that 

0"$(e) = -^(S^) (9) 

Then, the conditions of consistency of equations IpTjl 
and JSl can be written as: 



d x {ui*d v u u ) = d»{u x d u u v ) 



(10) 



In general case, there are 6 independent equations. 

Finally, the relativistic hydrodynamics of quasi-inertial 
flows is described by the equations © and H10J1 for the 
hydrodynamic velocities , and the equations (0 for the 



energy density: one should use derivative 



e+p{e) 



= *'(£) 



at any EoS p = p(e) to find function s(x). A serious 
problem is, however, to find non-trivial solutions for the 
field u^(x) of hydrodynamic 4- velocities. 



III. GRADIENT-LIKE VELOCITY ANSATZ 

One can try to satisfy to Eqs. i|5)l. (|lU[l for velocity pro- 
file by a use of gradient-like representation for it, namely, 



with condition of normalization, u^u^ = 1: 



(ii) 



(12) 



Then one can check that 10 is satisfied automatically, 
and H10|) leads to: 



d x (d' 1 (f>- □</>) = d' 1 (d x (f>- n<f>) 



(13) 



Thus, gradient-like velocity ansatz Hll|) reduce the prob- 
lem to equations (|T2"|> . (fHfl) . 

One can see that the above equation can be, in partic- 
ular, reduced to: 



Ucj> = F{<j>) 



(14) 



with any real function F that have to be solved together 
with (U2J. Note that if F(4>) = a + b<j) then ^ is the lin- 
ear inhomogeneous partial differential equations (PDE) 
and its any solution is a partial solution </>^ of inhomo- 
geneous PDE, plus general solution (f>h of correspondent 
homogenous PDE (if 6^0): 



and 



cj> ih ~ (t 2 -x 2 ),6 = 



d 4 P s(p 2 - b)f( P y 



(15) 



(16) 



where f(p) is arbitrary function with properties f*{k) = 
f(—k). Then the problem is reduced to a solution of the 



nonlinear integral equation i|12|) for f(p). If a = b = 0, 
the only potential <p = c+cgt+^CiXj (i = 1, 2, 3) with the 
constrain on the constants c„ : c§ — J2 C 1 = 1 is satisfied 
to these equations. It describes a relativistic motion of 
a medium as the whole. It is an open problem whether 
there are analytical solutions at a ^ and/or b ^ 0. 

The known quasi-inertial solutions correspond to 
F(<j)) = n/4> in Ea. H14|) . The value n = 1 generates gra- 
dient ansatz 4> = \Jt 2 — z 2 that gives the (1+1) boost- 
invariant Bjorken expansion along axis z, v = z/t; n = 2 
leads to cf) = yjt 2 — x 2 — y 2 , and, correspondingly, to the 
two-dimensional (1+2) Hubblc-like flow with cylindrical 
symmetry; at n = 3 one can get solution of l|14f> for cf> in 
the form d> = 



y 2 — z 2 describing spher- 



ically symmetric Hubble flow w M = x^ jr. The equation 
@ has the form d fl ^(e) = nx^/r 2 where number of space 
coordinates is equal to n. Then the energy density is de- 
scribed by the following expression 



e(r) 



de 



p(e) t 



(17) 



IV. RELATIVISTIC ELLIPSOIDAL SOLUTIONS 

One more possibility to satisfy to Eqs. or Q) be- 
sides of the gradient-like flows is to suppose a constant 
pressure in the EoS: p'(s) = const. This possibility was 
first used in as physically corresponding to a ther- 
modynamic state of the system in the softest point with 
the velocity of sound c 2 — 0. Such a state could be as- 
sociated with the first order phase transition. In A+A 
collisions it corresponds, probably, to transition between 
hadron and quark-gluon matter at SPS energies. The so- 
lution proposed in |8( has the cylindrical symmetry in the 
transverse plane and the longitudinal boost invariance: 



t x y z 
= y(-,v~,v-, -), 



(18) 



where r = \Jt 2 — z 2 , 7 = (1 — v 2 ) 1 / 2 and r is transverse 

(19) 



radius, r — \J x 2 + y 2 , 



1 + OLT 



describes axially symmetric transverse flow. 

The above solution has, however, a limited region of 
applicability since the boost invariance is not expected at 
SPS energies and can be used only in a small mid-rapidity 
interval [lOj , it is not reached even at RHIC energies 111. 
Most important, however, is that in non-central collisions 
there is no axial symmetry and, therefore, one needs in 
transversely asymmetric solutions to describe the elliptic 
flows in these collisions, e.g., coefficients. Now we 
propose a new class of analytic solutions of the relativistic 
hydrodynamics for 3D asymmetric flows. 
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First we construct the ansatz for normalized 4- velocity: 



{ 



a k (t)x k 



=} (20) 



where the Latin indexes denote spatial coordinates and a, 
are functions of time only. In this case a set of nonequal 
Qj induces 3D elliptic flow with velocities Vi = Oi(t)xi/t: 
at any time t the absolute value of velocity is constant, 
v 2 = const, at an ellipsoidal surface ^ a? x\ = const. 
Note that this solution is not gradient-like, so we follow 
in the specific way of a further analysis starting from J2J • 
The condition © of accelerationless in this case is re- 
duced to the ordinary differential equation (ODE) for the 
functions ai(t): 



da,i 
~~dt 



t 



the general solution of which is: 



cn(t) = 



t 



t + Ti 



(21) 



(22) 



where Tj is some set of 3 parameters (integration con- 
stants) having the dimension of time. The different val- 
ues Ti, Ti and T3 results in anisotropic 3D expansion 
with the ellipsoidal flows. 

The equation (JSJ is satisfied since we assume the con- 
stant pressure profile: p = Pa = const. The next step is 
to find solution of Eq. J3J for energy density e. Taking 

E &t/r, where 



into account that d^u^ 



(23) 



one can get 



(e+po)^^) +td t e + Y,a i (t)x i d i s = 0. (24) 



General solution of the equation is 



£ +Po = 



xi x 2 x 3 
t+Ti ' t+T 2 ' t+T 3 



) 



(t + TiXt + TaXt + Ta) 



(25) 



where F £ is an arbitrary function of its variables. If one 
fixes the parameters Ti that define the velocity profile, 
then the function F £ is completely determined by the ini- 
tial conditions for the enthalpy profile, say, at the initial 
timet = 0: e{t = 0, x) + p Q = F E (^ , §, %)/TiT 2 T 3 . 

If some value, associated with a quantum number or 
with particle number in a case of chemically frozen evo- 
lution are conserved then one should add the corre- 
sponding equation to the basic ones. Such an equation 
has the standard form [l2| : 



nd v u l 



u^d^n 







(26) 



where n is associated with density of the correspondent 
conserved value, e.g., with the baryon or particle densi- 
ties. A general structure of this equation is similar to 



what Eq. (|24|l has and, therefore, the solution looks like 
as (I2I3I1 : 



■<'i 



t+Ti ' t+T 2 ' t+T 3 



(i + T!)(t + T 2 )(< + T 3 ) 



(27) 



where the function F n is an arbitrary function of its ar- 
guments and can be fixed by the initial conditions for 
(particle) density n: n(t = 0,x)TiT 2 T 3 = F n (§h, §§, §§)• 

To establish a behavior of other thermodynamic values 
we use link between different thermodynamic potentials 
e = Ts—p+ un and utilize the thermodynamic equations 
based on the free energy density f(n,T) = e — Ts = 
fin — p. Since the volume is fixed (it is unit) the free 
energy depends on T and n only, df = —sdT + fidn, and 
the chemical potential fi = f \ n \T=const an d the entropy 
density s = -f, T \n=const- 

In a case of chemically equilibrated expansion of the ul- 
trarelativistic gas when the particle number is uncertain 
and is defined by the conditions and parameters of the 
thermodynamic equilibrium, e.g., by the temperature T, 
the chemical potential fi = (we suppose here that there 
is no other conserved values associated with charges, or 
the corresponding chemical potentials are zero or close 
to zero). Then / = — po = const, df = —sdT = that 
means the temperature T — const for such a system and 
the entropy s — (e(i,x) +po)/T where e(t,x) is defined 
by (23. 

If chemically frozen evolution takes place, the chem- 
ical potential associated with conserved particle num- 
ber is not zero and describes the deviation from chem- 
ical equilibrium in relativistic systems. The solution 
of differential equation nf n \ T=const — po = f(n,T) is 
f(n,T) = nc(T) — pa, where c(T) is some function of 
the temperature. Then it follows directly from the ther- 
modynamic identities that 



e{t, x) + po = n(t, x)(c(T) - Tc'{T)) 



(28) 



Since the structures of general solutions for n and e are 
found as l|25|) and (|27|l . the temperature profile has the 
form 



T(f,x) = F T ( 



xi 



3'2 



3'3 



t + Ti ' t + T 2 1 1 + Ti 



) 



(29) 



where Ft is some function of its arguments that is de- 
fined by the initial conditions for e and n and also by 
EoS £ = e(n,T). The latter can be fixed by a choice of 
the function c(T) in Eq. 1)28(1 . If the initial enthalpy den- 
sity profile is proportional to the particle density profile, 
K(f , f| , f ) ~ ^e(f , § , f ), then T = const (and so 
fi = c(T) = const) during the system's evolution for any 
function c(T) except the linear one: c(T) = a — bT when 
T is not defined by the equation l(28|l . In the last case 
n = (e + po)/a, s = bn. In another particular case which 
corresponds to EoS e+pa = anT with c(T) = — aTln(&T) 
one can get: 

T(t,x) = (e+p )/(an), s(t,x) = an(\n(bT) + 1) (30) 
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V. GENERALIZATION OF THE HUBBLE-LIKE 
FLOWS 

Let us describe some important particular solutions 
of the equations for relativistic ellipsoidal flows. If one 
defines the initial conditions on the hypersurface of con- 
stant time, say t = 0, then t is a natural parameter of 
the evolution. Such a representation of the solutions sim- 
ilar to the Bjorken and Hubble ones with velocity field 
Vi = a,iXi/t has property of an infinite velocity increase 
at x — > oo. A real fluid, therefore, can occupy only the 
space-time region where |v| < 1, or f 2 > 0. To guar- 
antee the energy-momentum conservation of the system 
during the evolution, all thermodynamic densities have 
to be zero at the boundary of the physical region, oth- 
erwise one should consider the boundary as the massive 
shell 0. Hence in the standard hydrodynamic apprach 
the enthalpy and particle density must be zero at the 
surface defined by |v(t,x)| = 1 at any time t. One of 
a simple form of such a solution (for the case of parti- 
cle number conservation) can be obtained from (|25|l . (|27() 

choosing F €l n ~ exp (- 



h 2 1 



e(t) +pa 



Iliit + Ti 



exp 



■ exp 



t 2 



<Y 2 



(31) 
(32) 



where f is defined by (|23|l . and the constants C e , C n , 
b e and b n are determined by the initial conditions as de- 
scribed in the previous section. As one can see, the en- 
thalpy density tends to zero when |x| becomes fairly large 
approaching the boundary surface defined by |v(t, x)| = 
1, in the other words, when f — ► 0. Thus the physically 
inconsistent situation when massive fluid elements move 
with the velocity of light at the surface f = is avoided. 
Of course, in such a solution one has to put a constant 
pressure to be zero, po = 0. 

As it follows from an analysis of a behavior of the ther- 
modynamic values in the previous section, the tempera- 
ture is constant if b e = b n , otherwise one can choose the 
temperature approaching zero at the system's boundary, 
e.g., for EoS which is linear in temperature, the latter 
has the form 



T 



const 



(33) 



0, \v(x)\->l b E >b r 



according to (pTO^l . 

Note that in the region of non-relativistic velocities, 

2 2 

v 2 = Y^, -fr 1 "C 1 the space distributions of the thermo- 
dynamical quantities (|31|) . l|32|) has the Gaussian profile: 



Po 



n t (*+ T .) e 



! 

, 2 2 ®i 



(34) 



The forms of solutions l|34|) are similar to what was found 
in Ref. ^| as the elliptic solutions of the non-relativistic 
hydrodynamics equations. In this sense the solution pro- 
posed could be considered as the generalization (at van- 
ishing pressure) of the corresponding non-relativistic so- 
lutions allowing one to describe relativistic expansion of 
the finite system into vacuum. 

One can note that the case of equal flow parameters 
Ti = and 6 £ = b n = induces formally Hubble-like 
velocity profile with the behavior of the density and en- 
thalpy similar to l|17l) at n=3, p — po, and with the sub- 
stitution t — ► t. 

The direct physical generalization of the Hubble so- 
lution for asymmetric case should be associated with 
the hypersurfaces of the pseudo-proper time f rather 
than with time t, that eliminates the problem of infi- 
nite velocities: v 2 — ^ < 1 at any hypersurface 
f 2 = const > 0. ft can be reached if one chooses the 
function F e and F n in (|25Jl . (|27|l in the form 



F 



(35) 



Then the generalized Hubble solution is 



{ 



t a\x l a-ix 2 a^x" 



e+p a = C t 



C 



T T 

aia 2 a 3 



h 



aia 2 a 3 



(36) 



where r = y/t 2 -'Eajx?, a, = a t {t) = t/(t + T t ) 
and constants are: C E = T 1 r 2 T 3 (e(0, 0) + p ), C n — 
7\T 2 T3n(0, 0). Again, the proportionality between e and 
n results in the temperature to be a constant during the 
evolution. If all parameters Ti are equal to each other, 
then Oj are also equal and solution (|3l)| just corresponds 
to spherically symmetric Hubble flow (at constant pres- 
sure) and t is the proper time of fluid element, r = r. 
Note that comparing to the standard representation of 
the Hubble solution the origin of a time scale is shifted, 
t — > t + Ti , and therefore the singularity at t — is ab- 
sent. Thus, if this solution is applied to a description 
of heavy ion collision, Ti should be interpreted as the 
initial proper time of thermalization and hydrodynamic 
expansion to which the origin of a time scale is shifted, 
typically t = Tj ~ lfm/c. As to a general case of an 
asymmetric expansion, the minimal parameter T t can be 
considered as the initial time t (at x = 0) of the be- 
ginning of the hydrodynamic evolution. In analogy with 
the Hubble flow the initial conditions in asymmetric case 
(|36|l can be ascribed to the hypersurface a : t = const 
so that tfj (x = 0) = 0. Note, that such a hypersur- 
face at |x| — ► oo tends to the hyperbolical hypersurface 
t = const since ta-(x) — > oo in this limit and so all a% — ► 1. 

The boost- invariant (1+1) solutions are also contained 
in general ellipsoidal solutions l|25|l , l|27|l for quasi- inertial 
flow. To get it one has to choose functions F e and F n 
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in the form (|35|l with another power: 3 — > 1; it leads 
to the same form of solution as (|36[) with replacement 
r 3 — > r . The next step is to suppress the transverse 
flow by setting T\ — > cxd, T 2 — ► oo (as usual, xi and 
CC2 denotes coordinates in the transverse plane and 23 
is the longitudinal axis), while the parameter T3 is fi- 
nite. This limit approach gives us a\ = a-i = and 
t — ► t = ^/(t + T3) 2 — x§ and results directly in the 
Bjorken solution at T3 = 0. Since Ti is a shift of a time 
scale to the beginning of hydrodynamic expansion, it is 
naturally to consider T 3 7^ as this was discussed above 
for the Hubble-like solution. This value transforms as 
T 3 — ► T3 = T3/7 at Lorentz boosts along axis £3. 

Note, that if one does not change the power 3 — > 1 
in H35fl and proceed to the limit directly in the equation 
(|36[1 . the particle density (and enthalpy) behavior will 
differ from the boost-invariant one as the following: 

— — ~* -('-ot)" < 37 > 

It is also the solution of (1+1) relativistic hydrodynam- 
ics at p — const but it obviously violates the boost- 
invariance: the particle and energy densities are not con- 
stants at any hypersurface and their analytic forms are 
changed in new coordinates after Lorentz boosts. 

It is worthily to emphasize that the physical solu- 
tions with non-zero constant pressure have a limited re- 
gion of applicability in time-like direction: if one wants 
to continue the solutions to asymptotically large times, 
then (e + po)t^oa ~ § — * 0, and this results in non- 
physical asymptotical behavior e — > —po, unless we set 
Po — 0. Therefore, it is naturally to utilize such kind of 
solutions in a region of the first order phase transition, 
characterized by the constant temperature and soft EoS, 
c s = dp/de w 0, or at the final stage of the evolution 
that always corresponds to the quasi-inertial flows. 

VI. CONCLUSIONS 

A general analysis of quasi-inertial flows in the rela- 
tivistic hydrodynamics is done. The known analytical so- 
lutions, like the Hubble and Bjorken ones, are reproduced 
from approach developed. A new class of analytic solu- 
tions for 3D relativistic expansion with anisotropic flows 
is found. The ellipsoidal generalization of the spherically 
symmetric Hubble flow is considered within this class. 
These solutions can also describe the relativistic expan- 
sion of the finite systems into vacuum. They can be uti- 
lized for a description of the matter evolution in central 
and non-central ultra-relativistic heavy ion collisions, es- 
pecially during deconfinement phase transition and the 
final stage of evolution of hadron systems. Also, the so- 
lutions can serve as a test for numerical codes describing 



3D asymmetric flows in the relativistic hydrodynamics. 
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